Dynamics of axialized laser-cooled ions in a Penning trap 
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We report the experimental characterization of axialization — a method of reducing the magnetron 
motion of a small number of ions stored in a Penning trap. This is an important step in the 
investigation of the suitability of Penning traps for quantum information processing. The magnetron 
motion was coupled to the laser-cooled modified cyclotron motion by the application of a near- 
resonant oscillating quadrupole potential (the "axialization drive"). Measurement of cooling rates 
of the radial motions of the ions showed an order-of-magnitude increase in the damping rate of 
the magnetron motion with the axialization drive applied. The experimental results are in good 
qualitative agreement with a recent theoretical study. In particular, a classical avoided crossing 
was observed in the motional frequencies as the axialization drive frequency was swept through the 
optimum value, proving that axialization is indeed a resonant effect. 

PACS numbers: 37.10.Rs, 03.67.Lx, 41.20.q 



I. INTRODUCTION 

The first proposal [1 for the use of trapped ions for 
quantum information processing (QIP) has been realized 
in its original and in other forms in several groups around 
the world [2[ [31 IH [5] . Recent work on large-scale entan- 
glement and scalability [6^, T of ion-trap quantum com- 
puting and the demonstration of the essential feature of 
quantum error correction [8 mean that trapped ions offer 
one of the most promising systems in which to investi- 
gate QIP in the short- to medium-term. All of the exper- 
imental demonstrations to date have been in RF traps of 
various designs. Although there are many advantages to 
using such traps, they have the disadvantages of the mo- 
tional heating of the ions due to their micromotion, and 
the relatively strong effects of patch potentials formed on 
the electrodes, particularly with smaller electrode dimen- 
sions. 

The use of charged particles in Penning traps for quan- 
tum information processing is being investigated by a 
number of groups [9l[l0l[ni[l2lll3l[l4l[l5l[l6]. By 

contrast with RF traps, there is no inherent heating of 
the ions by the trapping fields since only static fields are 
employed. The use of static fields also means that the 
trap can be designed to have larger trap electrodes, fur- 
ther away from the ions. Thus, the effect on the trapped 
particles of patch potentials on the surface of the elec- 
trodes should be smaller. A final advantage of Penning 
traps is the shielding from external magnetic-field fluc- 
tuations through the Meissner effect when a supercon- 
ducting magnet is employed; magnetic field stability has 
been found to be a limiting factor in some other experi- 
ments pTj . 

A major disadvantage of using ions in a Penning trap 
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is that the laser cooling process is greatly complicated 
by the large Zeeman splitting of the levels (much greater 
than a typical laser linewidth), and the overall negative 
energy of one of the two radial motions, the magnetron 
motion. This latter problem requires that, in order to 
"cool" the ions, energy must be added to the magnetron 
motion and simultaneously removed from the other two 
characteristic motions: the axial motion and the (radial) 
cyclotron motion. We have previously demonstrated a 
technique called axialization for overcoming this prob- 
lem using laser-cooled magnesium ions [18j. This is a 
variant on a technique previously used in conjunction 
with buffer-gas cooling in a Penning trap [19]. 

Although the hyperfine structure of ^^Mg+ could, in 
principle, be used for the storage of quantum informa- 
tion, the ground state and the metastable D5/2 level of 
^^Ca+ have already been shown to be a viable system 
for this purpose [20 . For this reason, and because of the 
availability of an all-solid-state laser system for address- 
ing the cooling and qubit transitions, we have chosen to 
investigate the laser cooling of ^^Ca+ ions. Doppler cool- 
ing of this ion in a Penning trap is reported in Ref. [9]. 
We present here the first demonstration of axialization 
of this ion in the presence of laser cooling, and the mea- 
surement of the damping rates of small numbers of ions, 
both axialized and unaxialized. 



II. ION MOTION IN THE PENNING TRAP 

The ideal Penning trap has the same geometry as a 
Paul trap, namely a pair of end-caps and a ring electrode, 
with all three electrodes having hyperbolic cross-sections. 
Traps with non-hyperbolic electrodes can give a reason- 
able approximation to a pure quadrupole potential and 
are often employed. A DC voltage applied between the 
end-caps and the ring provides confinement in the axial 
direction but repels the ion in the radial direction. Ra- 
dial confinement is provided by adding an axial magnetic 
field, forcing the ion into cyclotron-like loops. 



The motion in the axial direction is simple harmonic 
motion in a potential well. The motion in the radial 
plane is an epicycle, composed of a high-frequency cy- 
clotron motion and a lower frequency drift arising from 
the cross product of the electric and magnetic fields. The 
cyclotron motion is modified from that which would oc- 
cur in the presence of a magnetic field alone, due to the 
electric field. The high- and low-frequency motions are 
referred to as the modified cyclotron and magnetron mo- 
tions respectively. 

A detailed description of ion motion in a Penning trap 
may be found in [211 [22] • In summary, the potential from 
the hyperbolic electrodes is given by, 

= Uoi2z^ -x^- y^)/R^ , (1) 

where Uq is the DC voltage and = Tq -\- 2zq is a ge- 
ometrical factor depending on the trap radius (ro) and 
the end-cap separation (2zo). The magnetic field, B, is 
in the —z direction. The axial motion is independent of 
the magnetic field, with frequency, 

"^ = V^- 

The radial component of the ion motion may most nat- 
urally be examined in a frame rotating with the ions at 
half the true cyclotron frequency (i.e. at eB/2m) [23]. In 
this frame the force due to the magnetic field is exactly 
canceled by the Coriolis force. The two motions are then 
found to be circular orbits with opposite directions and 
the same frequency c^i, where, 

Lol=io^jA-Lol/2. (3) 

Superpositions of the two motions give rise to elliptical 
and (for equal amplitudes of the two motions) linear tra- 
jectories in the rotating frame. This is analogous to the 
way in which circular polarizations of light may be de- 
composed into linear polarizations and vice versa. 

Rotation in the same direction as that of the frame 
is found to correspond to modified cyclotron motion in 
the laboratory frame (at iJc/2-\-LJi) while rotation in the 
opposite direction corresponds to magnetron motion (at 



III. LASER COOLING IN A PENNING TRAP 

Although laser cooling of bound atoms was first 
demonstrated in a Penning trap [24 , it is rather more 
difficult than in a Paul trap. A relatively modest mag- 
netic field of 1 T can give rise to Zeeman splittings of 
many gigahertz (see Fig. [T]). The splittings can be so 
large that is is impossible to use a single, modulated, 
laser to address all transitions from a given level. The 
set-up and maintenance of many narrow-linewidth, drift- 
free laser systems is a significant complication. 
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FIG. 1: Simplified energy level scheme of °Ca^ at a magnetic 
field of 0.98 T. 




FIG. 2: The laser cooling efficiency of the radial motions is 
determined by both the spatial profile of the cooling laser 
(Gaussian TEMoo in the typical case) and the lineshape of 
the cooling transition (Lorentzian for a cold ion). The laser 
scattering rate depends on the instantaneous slopes of the 
tangents to these two curves and is varied by changing the 
spatial offset of the beam from trap center and by detuning 
the laser from resonance. 



As well as the large number of lasers required, the 
choice of detuning from resonance and the (radial) spa- 
tial offset of the cooling beam from the center of the trap 
is also important. This was investigated in an analytical 
model in Ref. [25 and further developed in Ref. [26 . 

To model the laser cooling process in a Penning trap, 
it is important to take account of both the lineshape of 
the cooling transition and also the spatial proffle of the 
cooling laser (Fig. [2|. The ion position in the laser's 
spatial proffle gives rise to a force varying approximately 
linearly with position along the axis perpendicular to the 
laser. 

Fa = -2amy . (4) 

The ion's velocity along the direction of the laser beam 
leads to a variation of the force on the ion as the ion is 
Doppler shifted into or out of resonance with the laser 
frequency, 

Fp = -2(3mx . (5) 

Adding these two forces to the equations of motion and 
making the assumptions that a <C cj^ and <C ojc, leads 
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to new solutions with frequencies f 

UJc' + I 

a - pUm 

+ ^ 7, 



(6) 
(7) 



That is, the motional frequencies remain unchanged 
by the cooling. The cooling rates (the imaginary parts) 
are determined by a and [3. It can be seen that efficient 
cooling of both motions is not possible due to conflict- 
ing requirements: positive (5 (i.e. a red-detuned laser) 
for the modified cyclotron motion and negative jS (blue- 
detuned laser) for the magnetron motion. A compromise 
may be found where both motions are cooled, but not 
efficiently [27 ^ 



(8) 



This corresponds to a red-detuned laser and a small spa- 
tial offset such that the laser is along the direction of the 
magnetron motion. 



IV. AXIALIZATION 

Since it is relatively easy to cool the modified cyclotron 
motion, it would be desirable to deterministically couple 
it to the magnetron motion. In this way the amplitudes 
of both motions could be reduced without direct damping 
of the magnetron motion. This is the essence of the axial- 
ization technique. We have described how the two modes 
of motion may be thought of as counter-propagating cir- 
cular motions in a rotating frame. In order to couple 
two such modes in a potential well, it is necessary to 
compress the well along one axis. This will result in an 
initially circular motion evolving into an elliptical motion 
which becomes more and more eccentric until it is linear. 
The motion will then become elliptical once again until 
a circular motion reappears, with equal amplitude but 
opposite direction to the original motion. 

The required compression of the potential well (in the 
rotating frame) is achieved by adding a static quadrupole 
field. To obtain such a quadrupole field in the rotating 
frame a rotating field must be applied in the laboratory 
frame. A quadrupole field oscillating at uOc can be decom- 
posed into two counter-rotating quadrupole fields each at 
c<;c/2. In the frame rotating uJc/2 , one of these rotates 
at ix)c and has no effect on the motion. The other is the 
desired static quadrupole. 

An oscillating azimuthal quadrupole field to couple the 
motions was first proposed by Dehmelt and Wineland [28^, 
[29j. Although the motion described is classical, their de- 
scription of the effect of the coupling field is in quan- 
tum mechanical terms. Such a drive has been applied in 
the presence of resistive cooling [28] and buffer-gas cool- 
ing pTOj where it is referred to as "sideband cooling" [33] 
or "radial centering". 



A. Quadrupole drive in the presence of laser 
cooling 

We have recently performed a detailed calculation ex- 
amining the effect of the quadrupole (axialization) drive 
when one of the radial motions is laser cooled [26 . In 
the rotating frame, the additional force on the ion due to 
the static component of the axialization field is given by 



-2rf 



(9) 
(10) 



where Vq is the amplitude of the axialization drive on 
the ring segments and vq is the inner radius of the ring 
electrode. We introduce a parameter, e, to describe the 
axialization strength. 
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and write v 
the ion as 



- iy. Thus, we can re- write the force on 
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The effect of laser cooling is described in the laboratory 
frame, with the co-ordinate, u\ 
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In the rotating frame moving at cOr^ with v = ue 
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where we have included the axialization term. Rotational 
symmetry is broken by the term and the solution is 
in the form of elliptical motion with uniform decay. This 
can be written as a superposition of two circular motions 
in opposite senses and with different amplitudes: 
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The complete solution is described in [2F. We define uja^ 
the axialization-drive frequency and A, half the detun- 
ing from the true cyclotron frequency, i.e. cOa = ^c-^ 2 A. 
The solutions to Equation [14] have motional frequencies 
similar to the unperturbed frequencies, c^i, with a fre- 
quency shift Sq. The damping rates are obtained from 
the imaginary part of the solution, 70. 



= cji + (^0 + ho 



(16) 



The solutions give. 



So = ±^\/ N ^ Vn^ ^ Z:^^ , (17) 
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FIG. 3: When {e/ouif > M^ an avoided crossine: of the 
motional frequencies (soUd hnes) is predicted, shown here near 
the magnetron frequency. The damping rates (dashed Unes) 
of the two radial motions are relatively high and insensitive 
to the axialization-drive frequency. The red dashed line is the 
damping rate corresponding to the upper frequency while the 
blue dashed line corresponds to the lower frequency. 



important experimentally, where (e/cJi)^ M^. In this 
case the axialization rate is faster than the rate of damp- 
ing due to laser cooling. Both the magnetron and mod- 
ified cyclotron motions are cooled and the sensitivity of 
the damping rates to the axialization-drive frequency is 
low compared to the other regimes (Fig. [3|. This means 
that very precise tuning of the axialization drive is un- 
necessary and that the damping rates are insensitive to 
magnetic field fluctuations. An interesting feature of this 
regime is that an avoided crossing of the motional fre- 
quencies measured in the laboratory frame is predicted. 



V. EXPERIMENTAL SETUP 

The setup has been described in some detail in Refs. O 
l3Q] . We perform laser cooling of small clouds of calcium 
ions in a Penning trap using diode lasers. An important 
feature of the system is that the ring electrode of the trap 
is divided into quadrants to allow the application of the 
axialization drive. 



A. Laser systems 
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The parameters M and N are given by 
M = {2a- f3u;c)/2u;i , 



N = A^ 
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(18) 



(19) 



(20) 



M is a measure of the strength of the laser cooling 
while N is dependent on both the laser cooling strength 
and the amplitude and detuning of the axialization drive. 

Equation [T7| implies that there are two frequency com- 
ponents in the rotating frame. Recall that in the labora- 
tory frame, the motional frequencies are = cj^ ± (cji + 
^o). There are thus four frequency components observed 
in the laboratory frame: 
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with 6o taking the two values given by Equation (17). 

We can vary M by changing the cooling laser beam 
position and detuning from resonance. M is also affected 
by the repumper laser detunings as well as the intensities 
of all the lasers. These parameters will also affect N. 
The other free parameters in Equation (20) are easily 



controlled by changing the axialization-drive amplitude 
and detuning. 

We investigated three regimes of axialization-drive am- 
plitude theoretically [26 . One of these is particularly 



A pair of laser diodes in commercial systems (Toptica 
DLIOO) provides light on the cooling transition (S1/2- 
Pi/2 5 397 nm) while home-built extended cavity diode 
lasers at 866 nm and 854 nm (D3/2-P1/2 and D5/2-P3/2 
respectively) close the cooling cycle. Repumping from 
the D5/2 level is necessary due to a small amount of popu- 
lation of the P3/2 level by the cooling lasers (ASE light at 
393 nm) and subsequent decay to the long-lived D state. 

Two lasers are used on the 397 nm transition on ac- 
count of the large ground-state Zeeman splitting of 
27 GHz, making a single, modulated source difficult to 
produce. The two lasers at 866 nm each have sidebands 
imposed on them via an RF source added to the diode 
current. The two pairs of sidebands then cover the four 
a transitions. A single 854-nm laser is found to be suf- 
ficient to repump the small amount of population found 
in the D5/2 level. 

All of the diode lasers are locked with side-of-fringe 
locks to low- finesse cavities. The cavities are constructed 
from Zerodur and are tunable via piezo transducers ar- 
ranged in a re-entrant design to minimize thermal drift. 
The cavities are kept in insulated containers which are 
temperature-stabilized. Roughly linear drifts of between 
15 and 75 MHz per hour have been measured. (The cool- 
ing transition natural linewidth is 23 MHz for compari- 
son.) 



B. Trap with segmented ring electrode 

The trap itself (Figure [4| was first described in [31 and 
is housed in a vacuum chamber kept at ultra-high vac- 
uum (10~^ mbar or lower). It approximates an ideal Pen- 
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ning trap but uses non-hyperbolic electrodes with coni- 
cal cross sections. The ring electrode is split into four 
quadrants to allow the application of azimuthal dipolar 
and quadrupolar potentials. The magnetic field for ion 
confinement is provided by an electromagnet with a max- 
imum field of approximately 0.98 T and a long-term field 
stability of better than one part in 10^. With an end- 
cap voltage, [/q, of 4.0 V, this gives cj^ = 27r x 141kHz, 
cj^ = 27r X 23.9 kHz, uoc' = 27r x 348 kHz. The true cy- 
clotron frequency, c^c, is 27r x 376 kHz. 

In the work described below, the axialization drive 
and the dipolar excitation signal are provided by 
high-resolution function generators, allowing the de- 
termination of relevant frequencies to 1 Hz or better. 
Fluorescence-photon arrival times were measured with 
a Photomultiplier Tube (PMT) connected to a Time-to- 
Amplitude Converter (TAG) and digitized by a Multi- 
Channel Buffer (MCB) in a computer. 



C. Measurement of damping rates 

Motional frequencies and damping rates were mea- 
sured with an established technique This involves 
applying an excitation field near-resonant with the mo- 
tion of interest and performing photon-RF correlation 
measurements. As the frequency of the excitation field is 
swept through the motional resonance, a tt phase change 
occurs in the relative phase of the excitation signal and 
the ionic motion. This is a feature of any forced damped 
harmonic oscillator and the width of the phase change 
yields the damping rate of the motion. Efficient excita- 
tion was achieved by choosing a suitable geometry (dipo- 
lar across the ring electrode in the case of radial modes). 

As the ion moves in the trap, it moves through the 
gradient of laser beam intensity. Thus, the probability of 
emission of a photon is determined by the position of the 
ion and is related to the phase of the ion's motion. In 
this way, the relative phase of the motion with respect to 



FIG. 5: Typical plot of the phase change as the excitation 
frequency is stepped across the motional resonance (in this 
case the magnetron motion). The fit is an arctan function 
whose width gives the damping rate. The motional frequency 
inferred from the fit shows that there is a shift related to the 
amplitude of the axialization drive. Inset Sample correlation 
data from the MCB. The fit is according to Equation |23| The 
modulation is due to the ion motion while the initial phase of 
the modulation depends on the excitation frequency relative 
to the motional frequency. 



the excitation field may be determined by examining the 
phase of the photon-RF correlation signal [32 . A typical 
correlation signal is plotted in the inset of Fig. [5] and the 
following equation is fitted, 

/(At) = exp(-aAt) [h + |c| sin(cjexcit At - (/>)] . (23) 

The underlying exponential decay is due to the normal 
statistics of delays between photon detection events — 
longer delays become progressively less likely, ojexcit is 
the excitation frequency, (j) is the most important fit 
parameter and is the phase difference between the drive 
and the motion of the ion. 

Equation [23] is fitted to the correlation signals obtained 
as the excitation frequency is stepped across the motional 
resonance. A plot of the phase differences against exci- 
tation frequency (Fig. [5| then gives two important pa- 
rameters: the resonant frequency of the motion and the 
width of the tt phase change (i.e. the damping rate). 



VI. RESULTS 

The magnetron damping rate was measured as a func- 
tion of the axialization drive amplitude. Measurements 
were made with small clouds of roughly 10 to 20 ions. 
Between measurements, the laser frequencies were re- 
adjusted to return the fluorescence rate to that recorded 
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FIG. 6: The magnetron damping rate is seen to increase with 
the axiaUzation drive ampUtude for a small drive. As the 
two radial motions are more strongly coupled, the otherwise 
poorly cooled magnetron motion is reduced in amplitude. For 
comparison, the damping rate in the unaxialized case was 
typically less than 5 Hz. 



at the beginning of the measurement. This reduced the 
changes in the damping rate due to laser frequency drifts. 
A series of damping rate measurements is shown in Fig.|6] 

The unaxialized magnetron damping rate was typically 
2-3 Hz. This is lower than expected from theory, in part 
due to the size of the ion cloud. With more than one ion, 
space charge effects increase the size of the cloud and it 
spends less time in the laser beam. 

The magnetron damping rate was seen to increase with 
increasing axialization-drive amplitude. The error bars 
are derived from the fit to the phase plot. It is difficult 
to infer the exact relationship between damping rate and 
axialization drive from the graph since the damping rate 
is also dependent on the laser frequencies. Drifts of the 
four laser-frequency reference cavities were measured to 
be in the region of tens of megahertz per hour. A typical 
damping rate measurement required twenty minutes of 
data acquisition. 

As well as a measured increase in the damping rate, 
the ion cloud's radius, dominated by the magnetron am- 
plitude, is seen to decrease. This effect was less marked 
than in previous experiments with a single magnesium 
ion [TSl. Again, this is due to space-charge effects that 
come into play when more than one ion is trapped. The 
aspect ratio of the ion cloud gives a clearer indication of 
the effect of the drive for larger numbers of ions (Fig- 
ure [t]). 



A. Phase measurements 
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FIG. 7: The aspect ratio of the cloud rapidly decreases as the 
axialization drive amplitude is increased, before leveling off 
above about 200 mV. Images a through / correspond to the 
first six point in the plot. The representative error bars are 
derived from the imaging system resolution. 
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FIG. 8: A phase change is observed at two frequencies close 
to the unmodified frequency of the magnetron motion. The 
separation and position are determined by the axialization 
drive amplitude and detuning from resonance with the true 
cyclotron frequency. The shape is in excellent agreement with 
theoretical predictions [26] . 



the frequencies of the resonances are shifted and addi- 
tional phase changes appear, so that there are observed 
phase changes at two frequencies near the magnetron fre- 
quency and two near the modified cyclotron frequency 
(Equation 17). 

It is expected that a phase change will occur at 
each of these "dressed states". An example of the two 



phase changes near the magnetron motional frequency 
are shown in Fig. |8] A similar pair could also be ob- 
served near the modified cyclotron frequency. 



In the unaxialized case there is no coupling between 
the two radial motions, and two tt phase changes are 
observed, at the magnetron and modified cyclotron fre- 
quencies respectively. When the coupling is introduced. 



B. Classical avoided crossing 

The frequencies of the phase changes are predicted by 
the analytical model summarized in Section P^V A[ taking 
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FIG. 9: An avoided crossing is observed when the axiaUzation 
drive frequency is scanned across the true cyclotron frequency. 
The fit parameters correspond to a true cyclotron frequency 
of 379.5 kHz, an unaxialised magnetron frequency of 23.9 kHz, 
and {e/uji) — 5.6kHz 

account of the laser cooling rates in a Penning trap and 
the strength of the axializing field. We operated in a 
regime where {e/uoi)'^ ^ M^, i.e. the coupling rate due 
to the axializing field was large in comparison to the laser 
cooling rate. 

In this case of a strong axializing field (200 mV), we 
expect to see an avoided crossing of the two measured 
frequencies as we scan the frequency of the axializing 
field across resonance with the true cyclotron frequency. 
We confirmed this by making measurements of the fre- 
quencies from plots such as those in Fig. [S] and plotting 
them against the axialization drive frequency (Figure [9|. 
From a fit to the plot, we can infer that the true cy- 
clotron frequency was 379.5kHz. The value for (e/a;i) 
was 5.6 ±0.2 kHz. 

From a simulation of the trap geometry in the 



SIMION [34] package, we estimate that the theoretical 
value for {e/uoi) should be 8.2 ± 0.8kHz for a drive of 
200 mV. The discrepancy between the measured and cal- 
culated values is in part due to the measurement being 
made on a cloud rather than a single ion. 

VII. CONCLUSIONS 

We have demonstrated the axialization of calcium ions. 
The use of RF excitation coupled with time-correlated 
photon detection allowed the measurement of the damp- 
ing rates and motional frequencies. The magnetron 
damping rate could be increased by an order of mag- 
nitude using this technique. We have also measured the 
predicted avoided crossing of the magnetron frequency. 

The results are in excellent qualitative agreement with 
theory. The error in damping rate measurements was 
relatively large due to laser- frequency drifts. Improve- 
ments to the laser systems should allow quantitative in- 
vestigation of the relationship between axialization drive 
amplitude and motional damping rates in the future. 

The axialization of laser-cooled calcium ions is an im- 
portant step in the investigation of the viability of a Pen- 
ning trap quantum computer. Future plans include the 
measurement of heating rates and spectroscopy of the 
Si/2 to D5/2 quadrupole transition. If the expected low 
heating rates are confirmed, early experiments investigat- 
ing arrays of Penning traps will allow for scalable quan- 
tum information processing with Penning traps [HI [13] . 
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